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After a brief discussion of the need for predicting the emergency 
ascent trajectory of a submersible and the means by which a mathematical 
model of an ascent can be derived, the second order, coupled equations of 
motion for a vehicle with varying mass and center of mass are derived. 

The equations of motion are then solved by a numerical step-wise 
procedure which is ameanable to programming on a digital computer. 

Ascent trajectories are calculated using data from model experi- 
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are all that are necessary to determine an ascent through an undisturbed 
medium, but are insufficient once the vehicle is disturbed from its vertical 
plane . 
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IliTRODU'J TICM 



An important operating node of deep submersible s is the ascending 
and descending node. Unlike large high speed military submarines which 
ascend and descend dynamically by use of a combination of speed and append- 
age deflection, deep subnersibles depend partially and sometimes entirely, 
as in the case of the 'iitldSTE , upon buoyant ascent and descent. Of the as- 
cending and descending operations, the most critical is during an emergency 
ascent in which most or all of the jettisonable weights are removed from the 
vehicle in order to obtain as rapid an ascent as possible. 

The majority of deep submersibles fall into a class of vehicles 
characterised by low speed and minimal appendage streamlining. They arc de- 
signed for underwater research and rescue missions which do not, except for 
the travel mode of the Deep Submergence Rescue Vehicle, require large 
speeds , They are equipped with many and diverse appendages for accomplish- 
ing their missions such as manipulator arms, lights , propulsion motors, 
sampling equipment, TV cameras and mating skirts, all of which effect, to 
varying degrees, the streamlining and hydrodynamics of these vehicles. 

In order to study the motions and ascent trajectories of cub- 
mersibles in normal or emergency ascent without model or full scale tests, 
it is necessary to develop a mathematical model of the vehicle which can 
accomodate the effects of varying mass and center of mass. 

This work is an attempt to develop the mathematical model using 
a set of second order equations of motion, and solve these equations by 
use of a high speed digital computer. 

The primary sources of material for this work are the works of 
Abkowita, Dog an, Gcrtlcr and Hagen, Lamb and Strumpf, but the secondary 
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sources, acknowledged or not, axe no less important ot the development of 
this work. 

Tlie problem ox determining the notions of a totally submerged 
body has been treated by many authors, however, there appear to be only 
two basic means for arriving at these equations. These arc by energy 
methods and by vector calculus. 

The first method is that vised by Sir Horace Lamb in his book 
"Hydrodynamics" which first appeared in 1879 (see ref 1). The cnorgy 
method is based upon the existence of a single valued velocity-potential 
which implies that the motion of tho fluid is entirely due to that of the 
submersible, and is therefore irrotational and acyclic, This method leads 
to a completely general set of equations representing the rigid body dynam 
ics, but the hydrodynamic forces that result represent only tho so called 
added mass or added inertia effects and not all the forces acting on the 
body. Of tho forces that are missing the most important are those due to 
circulation, separation and vortex shedding. 

The second method of developing the equations of motion is to 
combine a vector expansion of Newton’s laws of motion, expressed as follow 



F = (momentum) - the vector force acting on the body 



( 1 ) 



j-i = 



dt 



(Angular momentum) - the vector moment acting on the 

body 



( 2 ) 



with a Taylor series expansion of tho hydrodynamic forces and moments ex- 
pressed as: 



Forces 

Foments 



f (Properties of body; Properties of motion, 
Properties of fluid) 



( 3 ) 



2 



which reduces to 



-L. = f (Properties of notion) 
M 



(4) 



for a particular vehicle in a particular fluid* Tills method has been used 
by A. Strumpf in i960 in developing Ms "Equations of Potion of a Sub- 
merged Body with Varying Mass" (see ref 2) and by Prof, M* A, Abkowitz in 
1949 in his lecture course and notes on "The Dynamical Stability of Sub- 
marines" (see ref's 3 and 10), It too results in a completely general set 
of equations describing the rigid body dynamics, but its greatest asset 
lies in the generality of the hydrodynamic forces and moments, The hydro- 
dynamic forces and moments obtained from a Taylor series expansion in- 
clude not only all the added mass effects but also the circulation, sepa- 
ration and vortex shedding effects. The only limitation on this method 
is the availability of theoretical or experimental data to use in the equa- 
tions. 

From the proceeding discussion and equations (l through 4), it is 
obvious that the derivation of the rigid body motion may be completely sep- 
arated from the problem of developing a suitable form for the hydrodynatiic 
effects. Therefore, the two methods of development may be broken down and 
part of each used. 

This work shall attempt to take advantage of the simplicity of 
the energy development of the rigid body dynamics, while retaining the gen- 
erality provided by the Taylor series expansion of tho hydrodynamic forces 
and moments, 

The solution of the equations of motion that will bo developed 
employs a technique suggested in 1964 by G , Parissis in his solution of 
linearized equations of motion for Leave and pitch of a ship (see ref 4) 
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end again in 1 96? by L. h, iicCloskcy for tho digital siniulation of the DSRV 
control system and autopilot (see ref 5)» 

The coefficients used to tost the solution of the equations of 



motion come from a series of model experiments conducted on tho DSRV, The 
results of the computer solution mill not, ho:: ever , be the ascent trajectory 
for the DSRV since there is an assumption, in this uork, of zero propulsion 
forces in the development of the equations of motion < 

The assumptions and the statement of the problem are given in Chap- 
ter I, the equations of motion are derived in Chapter II and a method for 
solving the equations of motion is developed in Chapter III. 
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CHAPTER I 



F0IMJIATIG1I OF THE l’KOBLIfilu 



The problem Involves a deep submersible vehicle of given geometry 
and dimensions which is in free ascent through a stationary fluid. 

1-1 Assumptions 

Tho vehicle is assumed to be a rigid body with no clastic defor- 
mations of a vibratory nature. It has six degrees of freedom, three trans- 
lational and three rotational. Ue shall bo interested in motions in both 
the horizontal and vertical planes of motion. Velocities aro small. Hydro- 
dynamic effects of second order in acceleration shall bo considered negligi- 
ble, Tho only body symmetry is port and starboard, 

1 - 1.1 

Tho neglecting of vibratory motions is reasonable since the fre- 
quencies of vibration of tho hull acti'ig as an elastic body aro of different 
orders of magnitude than the frequencies of motion and do not excito the 
latter. Elastic deformations due to vehicle compress ability must be in- 
cluded since they directly affect the buoyancy of the vehicle. 

I - 1.2 

Tho interest in both pianos of motion is duo to the desire for 
as general a sot of equations as possible. 



1-1.3 



The assumption of small velocities 



is realistic in that free as- 



cent velocities 



and open ocean currents aro in general of order less than 



ton knots , 
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' I - 1 .4 



The use of a sot of second order equations of notion is deemed 
necessary to appropriately describe the hydrodynamic cross-coupling that 
takes place in a problem such as this* The second order acceleration ef- 
fects are assumed to bo aero on the basis of potential theory (seo rof i), 
This, however, has not been experimentally verified 0 



1-1.5 



Neglecting as symmetries due to relatively small appendages which 



are not control surfaces, 



there are few vehicles which operate in the ocean 



environment which do not possess port and 



starboard symmetry. 



For this rea- 



son the assumption of port and starboard symmetry has been made „ 

1-2 Initial and Other Co n ditions 

The vehicle is to be initially at rest relative to the inertial, 
earth fixed, sores. The control effectors, propellers, thrusters, rudders, 
dive planes, etc., are inoperable and/or in neutral position. The driving 
force for the vehicle shall be a decrease in weight due to jettisoning of 
ballast or an increase in volume due to vehicle decompression as it rises. 

These conditions, though arbitrary, serve to restrict the problem 
to one of manageable size. 
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CHAPTER II 



EQUATIONS OP ACTION 



II - 1 General ; 

The derivation of the equations of motion for a rigid body in six 
degrees of freedom with a varying mss and center of mass follows closely 
that of Dr. Pierre Dogan in reference (5) t which is based primarily upon 
the development of the equations of motion by Lanb in reference (i). The 
derivation of the hydrodynamic force equations follows that of Professor 
Martin Ablcowita in reference ( 3 ) except that, in this work , the second 



ordor terms are retained. The forces due to gravity are determined using 
the transformations set forth in reference (6) and the angular velocity 
transformations given by Professor Abkouits in reference ( 7 ). 



The axis 



systems necessary to describe the notions 



of a body and 



its 



trajectory tlrrough a fluid include 



a body fixed system and an inertial 



system. 



II - 2 Axis System s 

The equations of motion for an ascending vehicle must be written 
in an earth fixed axis system in order to deter-inine the trajectory of the 
vehicle relative to some fixed point. An additional body fixed system is 
required in ordor to describe the hydrodynamic interactions between the 
vehicle and the water* 



II - 2,1 Inertial Ayj s System 

The right-handed earth fixed axis 
is assumed to bo an inertial system for the 



system, Ag (see figure 11,1 ), 
lee.son that the accelerations 
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Fi'-'-nro 1 C core! :i)i ate Syrtrns 



6~A 



of a point on tho surface of tho oarth arc an order of magnitude smallor 



than those vrliich aro of importance to tho motions of tho vehicle , An, is an 
orthogonal set of axes fixed relative to the earth such that components x^ 
and yp aro in a horizontal plane, and the axis is vertical and directed 
downwai’ds , 

II ~ 2,2 Body Axis System 

The right-handed body axis system, A^ (see figure II .1), is fixed 
in tho vehicle such that advantage is taken of the assumed principal plane 
of symmetry by placing tho origin of the system in this piano. The axes of 
this system ares 

x axis - the longitudinal axis, directed from tho after to tho 
forward end of the vehicle, 
y axis - the transverse axis, directed to starboard, 
z axis - tho normal axis, directed from top to bottom (deck to 

keel) . 



Tho xz plane is the assumed principal plane of symmetry, 

II “ 2.3 Body Axes Orientation 

Angular displacements of relative to A^ aro specified by a set 
of modified 'Euler angles' which aro taken as positive in the sense of 
rotation of a right-handed screw advancing in the positive direction of the 
axis of rotation. 

The orientation of Ap relative to A is described in terms of a 
roll angle $5, a pitch angle 0 and a yaw angle "f , Beforo defining these 
angles, an order of rotations must bo chosen since finite rotations are not 
true vector quantities and do not obey the rules for adding lectors. Tho 
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order chosen here conforms to that set forth in reference (6) viiich is: 



(1) rotate about the initial z - z r , axis tlirough an angle of 
yaw f , 

(2) rotate about the new position of the y y^ axis through an 
angle of pitch 0, 

( 3 ) finally rotate about the new position of the x ~ x axis 
through an angle of roll $5* 

In accordance with the order of rotations above wo have the 
following definitions : 

G - the angle of pitch; the angle of elevation of the x sods; th 
angle between the x axis and the horizontal plane o' ( -,y^, , 

"P - the angle of yaw; the angle fron the vertical plane to 

the vertical plane xz„, 

£j 

<f> - the angle of roll; the angle from the vertical xz P Pi ane to 
the principal plane of symmetry xz e 



The successive rotations required to specify the orientation of 

the body axes relative to the earth fixed axes can be described by throe 

orthogonal matrices [yj , I C , ml: 

L Jz- V -- J y 1 ' — 'X 




sin ~jr 
cos y 

0 



0 

0 

1 



cos G 0 ~ sin 0 




sin C 0 cos C 
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ML •• 



1 

0 

0 



0 

con (■) 
-sin <!> 



san p 



cos p 



Unit vectors in A n and A,, are related by an orthogonal matrix? 

IS Ji 



_i B 




3 t, 


^B 


- T 

B 




ks. 




k, 


B 







where is the product ox the three orthogonal matrices defining the 
rotations : 



B 



cos G cos f 

-cos sin} ,/? 

■:-sin p sin G cos V r 

sin <!> sin ~f 
+sin G cos ^ cos)'" 



cos 0 sin f 

cos $ cos Y' 

•:-sin p sin G sin'/' - 

-sin cos Y'" 

•!cos sin G sin 



-sin G 
sin cos G 

cos cos 0 



Velocities in and A^ are also related by the transformation 



matrix T^: 



u 




'k‘ 

X, 


V 


!! 


e 

y E 






* 


V 




. Z E_ 



The vehicle angular velocities in A n and A,., are, however, related 
“ B Ji 

by a non- orthogonal matrix, which is the sum of tliree components along the 
y and x axes of magnitude ^v. 6 and p , 
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■ 1 
p 




0 




0 




t; 

6 


q 


i 


0 


+ Mx 


• 

0 




f 

0 


r 

_ 


r 

r 




0 




0 



i 


0 


-sin 


e 






£ 


0 


cos $ 


cos 


G 


sin <!> 


c 


f 

G 


0 


-sin ^ 


cos 


0 


cos ^ 




r 

t 



T. 



* 

6 



f 



II ~ 3 C onservation of Momentum Equation s 

The derivation of the dynamical equations for a vehicle with a 
varying mass and center of mass have been treated in tyo different manners 



by: 



(1) Albert Strumpf of Davidson Dab, (see ref 2.) using vector 

» 

calculus, and including variations in mass, moments of in- 
ertia and CG position. 

(2) Pierre Dog an of HIT Instrumentation Lab, (see ref 5) t using 
a Langrangian formalism, including time variation of the 
inertial tensor, and making assumptions as to the form of the 



movable veiphts. 



The treatment by Strumpf assumes a body that is a rocket, so that mass is 
discharged from tlis body. The treatment by Dogan, for the Deep Submergence 
Rescue Vehicle (DSRV), assumes that the mass of the body is constant but 



allows the position of the center of gravity to change. 

The development by Dogan avoids the long vector..! cerjv..tic He- 



ll 



involved in Sti'umpf's development, irlrilo sacrificing generality by not 
including a change in mass. This feature could be included if it- irere 
desired and so the approach by Dog an v T ill be used. 

Nino generalized velocities and coordinates are sufficient to 
describe the motions of the vehicle. These are: six velocities (u, v, it, 

p, o, r) to describe the vehicle and three coordinates (x , y. , z ) to 
describe the motion of the center of gravity. Tho following six Lagrange 
equations will give tho needed vehicle momentum equations • (See Lamb, 
,, lIydrodynand:cs ,L 6th Ed, page 168) 
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v 



K 









r 9P ~ H 

•3T 
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* - N 



( 2.3 “ 
( 2.3 - 
( 2.3 - 
( 2.3 " 
( 2.3 - 
( 2.3 - 



1) 

2 ) 

3 ) 

h) 

5 ) 

6 ) 



T- is the vehicle total kinetic energy. X, Y, Z, K, H, N arb 
tho generalized forces and moments of which some can be further described 
by a potential function and others represent friction and drag. 

The three Lagrango equations necessary to describe the momentum 
balance equation for the jettisonable ballast subsystem are assumed to 
reduce to quasi- steady equations defining the Cl location from tho integral; 
of the various ballast release rates. 

This proceeduro is as much a necessity as it is a simplification 
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since, in tho cc.se of dropping ballast, tho ballast way take almost any 
form from liquid mercury and iron shot to large blocks of metal or pieces 
of equipment* The dropping of a liquid or a granular solid can be reason- 
ably modeled as a function of time but the dropping of chunks of metal or 
pieces of equipment would create singularities in a functional relation- 
ship* The obvious answer would bo to use a combination of a smooth func- 
tion and steps to obrain a reasonably accurate model of the deballasting of 
a vehicle * 



There is one additional factor which also affects the decision to 
reduce the function to a quasi- steady process* The total bail-last dropped 
is no more than three percent of the total vehicle weight and the rate at 



which it is removed is of order *3 percent of the total vehicle weight per 
second* This would then say that any Taylor series expansion of this 
function, which retained terms commesurate with the second order expansion 
to be used in obtaining the hydrodynamic forces, would contain at most the 
linear terms* 

i 

The effect of such an approximation is entirely dependent upon 



tho length of the time interval over which tho process is assumed to be 
steady and shall be discussed in conjunction with the computer program* 
Suffice it to say here, that until accurate model tests can be conducted, 
the effect of this assumption is truly unknown, but appears to bo of incon- 
sequential magnitude . 

The equations describing tho variable position of the CG in the 
body axis system are-* 
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(2.3 - 7 ) 




(2.3 - 8) 



(2.3 - 9) 



whore x , y , z arc the components of the CG of the vehicle with all of 
GB GB GB 

the ballast, W is the vehicle weight including all jettisonable ballast, 

x. , y. , z are the CG's of the N jettisonablo ballast weights h T . . 

-*• x i 3- 

Tho vehicle kinetic energy is the sura: 



T = | Vj •!• i-oI VJ (2.3 - 10) 

v ^ G u 

where n, V , and oe> are the vehicle mass, the C!G velocity, the inertial 
G G 

tensor about the CG and the angular velocity vector. Defining axes x* , y' , 
z' tlirough the CG parallel to the vehicle body axes, one has 



x'x* 


T 

X'y* 


i 

x'z' 






'y'x* 


w 


I 

y f z ! 




(2.3 •• U) 


Z * X* 


11 f.'y' 


I 






V'y' 


» « t c c have 


the usual meanings: 
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x'x' 



( y 



-r z'^) dll 



Y 



w - sss (xV,d 

V 



cm 



Computing the components of V :ln body axes one has* 

u 



% zz (u Q z g " ry G '* S G» v * rx G “ ps G * y G’ w •• KTq “ qX G * t } 

( 2.3 - 12 ) 

Substituting (2.3 ~ 11) and (2. 3 ~ 12) into (2.3 “ 10) one get 
the kinetic energy T and its partial derivatives: 






^ ~ = M (u qz ~ ry •!- x, ) 
3U 1 G b (J 



7>T r 

Tv" = B <v + ”*} - pz G V 



2> T 

A 

■3W 

VI 



- “ m (vr py - qz r * :',.) 



G 



? - J w p + :c X ' K ,r 



■JP x»x 



x y 



(2.3 - 13) 
(2.3 - 14) 
(2.3 - 15) 



nz (v rx. - pz •: . y„. ) 
u u Ct U 



+ ny G <W • > Py G " * %> 

B T v 

"I t t q •!' I , ,r + I p 

•^q y y y y'x ' 1 

- ^ (,: + p y r , - « x „ + 1 ) 



(2.3 - 16) 



+ TO G <u + qa G " ,,y G * V 



9T 



r~~ - X r I ,p I , ,q 

hr z*z* z'x’ z y : 

" r ' iy G (u ■' ( i'% ~ ry G ’• V 



(2.3 - 1?) 



+ “ti (V '• «b - P“g '• V 



(2.3 ~ 16) 
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Tho inertial tensor 1„ about tho CG can bo algebraically related to the 

G 

inertial tensor I about tho origin (defined in y, c) and the distance 

o 



V v *■ 



, . I 'or example: 

G 



1 - I , , in ( E ? y?) 

xx x'x‘ G ‘G 



I 1 
xy x'y’ 



’Vg 



Substituting theso for the inertial tensor , 1^, equations (2,3 
through ( 2,3 ~ 18 ) become: 

2 T v 

— — = I _p *i- 1 q •!- 1 r 
<>P xx 1 xy 1 xz 



- 16 ) 



t:iz g (v V 



•i- my (\r t •/. ) 
G u 



"2> T. 



—I = x q I r I p 

yz pc 



( 2,3 - 19 ) 



- mx (w 1- z,.) 

Li It 



mz_ (u x , ) 

U G 

>T 

— -- ~ I r + I p + I a 
za zx 



( 2.3 ~ 20 ) 



” y G U * V 



+ mx (v •!- y ) 
G G 



( 2,3 - 21 ) 



Substituting these equations and equations (2,3 ~ 13) through ( 2,3 - 15) 
(without primed subscripts) into the Lagrangian equations of motion 
(2,3 " 1) through (2,3 - 6), one gets the dynamical equations: 
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X - m [u - rv + qw - x (q 2 * r 2 ) + y G (pq - r) Z G (P r 

•!- 2qz - 2ry •!• xl (2,3 ” 22) 

u b bj 

Y = m [v - pw ru - y (r 2 + p 2 ) * Zq (qr - p) + * G (qP 

2rx G - 2pz c yj (2^3 “ 2 3) 

Z " M [tt - qu + pv - a G (p 2 ■!- q 2 ) •!- x Q (rp - q) + y G (rq p) 

2py G - 2q> G z G ] _ (2,3 •* 24) 

K - I Jp v I (q - pr) I (r * pq) + I (q 2 - r 2 ) * (I - I ) qr 

x>: xy 1 ^ xz vz ^ - V - Y 



xz 



yz 



~ mz G (v ~ pr ru + rx s- y^) 



*• wy (w - qu + pv - qx^ + 



-!• I p •!- I a I r + roc,. (Jr + y a ) 
XX xy * xe G G G 



M = X q •!- I (r - qp) I (p + qr) + I (r 
yy yz A yx zx 



(2,3 - 25) 



p ) * (I - I„„) rp 

XX ZE 



MX, 



ti 



(w - qu pv py ” ) 



•< nz G (u “ rv qu - ry G *5- ) 

c * • 

*!- I q '!• I r -1- I p •!- wy r (ip •!■ z_r) 

yy yz yx r v g 

“ 4/ * 4.x tp - **> * 4y w - + rp) * V (p ‘ 



2 



( 2,3 - 26 ) 



« 2) (I yy - I xx ) 1X1 



lay (u - rv -i qw •!- qz + x ) 
b b b 



hXq (v - pr ru - pL + y ) 



+ x 



r * I zx P '• r y; %<i + X^p) 



G G 



(2,3 - 27) 



1 ? 



where I , I , arc about axes through the vehicle origin. These 
xx xy 

* * e 

equations also contain derivatives of tho form 1 .1 . . . . . x x , . . * . » 

1 xx — >' 



*y 



U* G 



The tensor I can be thought of as being made up of two parts; a 

o 

constant part representing tho vehicle with all jettisonable ballast re- 
moved, X , and tho j o tti sonablo ballast, I, , such that: 

X *■* 



1 I •; I 0 

o 1 2 



( 2,3 - 28 ) 



The contribution of the ballast can be modeled by assuming that each weight 
that is part of the jettisohable ballast is lumped at its center of gravity 

(x . i y . . 55 . ) . 

£1 gx EX 

Then: 



I = *S m (y 2 . z 2 .) 

xx, 2 ^ x gx go. 



( 2.3 - 29 ) 



N 



1 ~ m.x .y 

xy, 2 i = i a- g:>- gi 



(2.3 - 30) 



N_ 

r ,= m. x . z . . . . 

xx, 2 i'S"! 3 - g ;L g'J- 



( 2,3 - 31 ) 



where II is the number of weights. 

1 is the time varying part of 1^, however, to bo consistent with 
the quasi- steady approximation made in developing equations (2.3 ■■ 7 
through 2,3 ~ 9) i a quasi-steady change in the inertial tensor must also 
be assumed. This then .says that tho terns involving time derivatives of 
the inertial tensor and the center of gravity can be dropped, since they 
are aero during the time interval ever which tho process is assumed to be 
steady. 



The equations can be further reduced whon it is recognise 
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duo to the assumption of an xz plane of symmetry and a further assumption 
that the ballast will bo dropped cyasac trie ally, tho products of inertia. 

I « I , « I „ X « X = I » 1 r.- I :: 0 

xy.J- yx.l yz,l zy,l xy,2 yx,2 yz»2 zy,2 

Thus, tho total moment of inertia becomes quasi*- steady and can bo used in 
the equations of motion without 0, i or 2 subscripts o The final form of 
the dynamical equations then becomes £ 

X -- m £~u - rv -!• qu " x^ (q 2 •!- r 2 ) y G Cpq - r) z G (pr '• q)j 

(2.3 “ 32) 

X — m [v - pi-r *:■ ru - y G (r 2 d p 2 ) d (qr - p) ; ■ (qp d r)J 

(2.3 - 33) 



= m 0 



xt - qu pv 



z G (p 2 t q 2 ) x Q (rp 



q) y G (rq + p)j 



(2.3 ** 3 / t) 



K ~ I £ d 1 (r pq) (X -I ) or 
xx xs zz yy x 



d m Q- s G (v - p-.r v ru) d yn (u d pv - qu)) (2.3 - 35) 



H - I q X (r 2 - p 2 ) •: (I - I ) rp 
yy sx xx 2s 

d m f- x G (it - qu d pv) 2 .^ (u d qu - rv)J 



(2.3 - 36) 



N = 1„ r d I (p ~ rq) (I - I ) pq 
zz ex ' yy xx 

in H-y^, (u - rv d qu) d x G (v ru - pu/j (2,3 - 3?) 

where n, x , y , s , X , } , X , X are all quasi- steady functions of 

u It Li XX yy 7»Z ^ /i 

time , 

The forcing forms for tho hydrodynamic equations are made up of 



19 



gravity forces, hydrodynamic forces and. propulsion induced forces. Sine© 
tlie problem has been defined as a free ascent problem thoro irill be no 
propulsion forces included. Because of this restriction, the experimental 
hydrodynaraic coefficients should bo obtained without propellers running. 
The total forcing terms are then: 



X - X G ■: X H 

1 - h + Y l, 

Z = Z * Z H 



K ss K g K 



H 



11 “ I: G + \ 

K = H q v K h 



(2.3 38) 



IX - 4 Gravity F o rces 

The hydrostatic pressure field induced by gravity creates- a 
buoyancy force B through the OB. This force varies with ambient water 
density and vcliicle volume. The instantaneous weight cf the body, VI, act- 
ing through tho CG, is made up of VI , the weight of the body without jetti- 

ji 

sonable ballast, less W, , the weight of the ballast components removed. 

3- l]_ 

During the initial phase of the ascent the items in V4 

i = 1 

are increased until all of the ballast components have been removed, at 
this point tho buoyant force (-B •!- V/) becomes a maximum, if we neglect 
changes in vehicle volume and density. This maximum force is sustained for 
the rest of the ascent. 

The instantaneous buoyant forco can thus be represented by: 

M 



b + vj - s v;. 

i^i x 



(2.4 - 1 j 

\;hcro K ~ the number of ballast components removed. This force rots 
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upwards along tho local vertical. Dae to tho choice of origin of the body 
axes, tho gravity induced torque is "r r x. B whore "r is the CG position 
vector in body axes raid B is the vehicle buoyancy vector acting up along 
the local vertical. Resolving along body axes tho gravity forces become: 



X n --- - (v/ - B) sir. 0 
= (VJ - B) cos 0 sin {> 
(VI - B) cos G cos /> 






Z r 



(2,4 - 2) 
(2.4 - 3) 
(2.4 - 4) 



Fy, i= (y_,W - yj ,B ) cos 0 cos < r ' •• - z^B) cos 0 sin {> (2.4 - 5) 



B 



N 



(: • V7 •• x^B) cos 0 cos </> - (eqW - s^B) sin 0 (2.4 - 6) 



N -- (xdi - x-sB) cos 0 sin f5 (yftf - y p B) sin 0 



G 



(2.4 " 7 ) 



11-3 Hydrodynamic Forcer. 

Tho hydrodynamic forces and moments that act on a body moving 
through a real fluid are the result of : 

(1) hydrodynamic inertial affects (linear added mass terms), 

(2) skin friction, separation arid cross-flow drag effects, 

(3) circulation offects. 



Those offects are all functions of the velocities and accelerations of tho 
body. Therefore, the hydrodynamic forces and moments can bo expressed 
functional) II y o.s • 



I'll 



I 



H 



” f («. V, V, P, q, r, u, v, w, P, q, v) 



( 2.5 - 1 ) 



This function nay bo rcducod to a workable form by expanding tho 
function in a Taylor series. Expanding the function in this form 1 squires 
that the function yrd its derivatives bo continuous in the. region of the 



values of the variables under consideration. A topical second order expan- 
sion of one of the force equations would be of the forms 



Y_v •- „ <> x v 

A ~ A V u — — * v • co« 6 oeroeace»ooeoc { J ^ - ' 

o Y*u i,v 



, . *X . .. *X 

v u ’» v r " r r y * t o c o 

Mi 






e o c o e c * e t fc c o 



-i T 






, , <„2 A ... 2 & .. 

• 2 tU n- • v • « « 

_ 1 >U \m 



L 2 b 2 X 

_JCtO€‘Ot»C » i. ’ 

'V 

2 ” 2 
t X ^ X 

'i' UV — — 'r U 1 T I ccccceeette U. 1 * 

2W-J c?u>r 



*dr 

^x 



* 2 X , * 2 X ^ * 2 X 

r pq — v Pty^ + v r<i»r 



. >2x . /'X , u « A 

*!- UU W + UV t « < C . C » t ur 

^)U2|U SU2>\/ 



• 


5) A _ 


*r ril 


r & 




2 . 


• c 


b X 


+ uu 


— tz + 




2>u 



. -& 2 X 
-*r — ~ r 
2>rar 



etcccccccccc 1 



rr 
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£ 



, .. a 2 X , .. * 2 X 

*** pq • ' •’ cccc.eccec.e.c.oco,,, t qr 1 ~G* 

>p*q ^qfr 



(2.5 - 2) 



This equation contains ninety-one constant, linear and socond order terms 
arising from the Taylor series expansion. The number of terms, however, 
can be reduced to thirty-throe on the basis of the problem assumptions and 



restrictions delineated in Chapter T , 



The constant term ts dropped to conform with the requirement that 



initially tho only disturbance is duo to a buoyant force xdrich is included 
as a gravity tern, 

Tho terns involving products of accelerations viith velocities or 
accelerations are dropped because the second order effects vroro restricted 
to volocitios only. In addition the results of potential theory indicate 
that these derivatives are aero, (See ref 2) 

After eliminating all but the linear terms and the second order 
velocity terms, the Taylor series expansions of the hydrodynamic forces and 
moments become of the forms 



X„ = X U X V * X XT X p X q X r 
H u ” ,r " - 1 



vr 



+ X.u + X„y X.W X.p •!• X.q X.r 

u v vr • P q v 

l V 2 I -V' 2 , y 2 , y 2 , V" 2 , y 2 

V X 2 u A 2 v V X 2 u d- X 2 P V X 2 q X 2 r 



•!- X uv d X uxr X up d- X uq d X ur 
uv ux-r up uo ur 



d- X vxr d- X vp X vq d- X 'vr 
VX’7 vp vq vr 



d* X vp X xrq d- X v:r 
vrp ^ xrq 1 vr 



•!- X pq X pr d* X or 
tl pr qr 



(2.5 - 3) 



The terms in this expression are seen to fall into ono of three 
catagorios, namely: 



(j ) added mass or inertial, 

(2) second order non- inertial, 

( 3 ) linear. 

In order to further reduce the number of terras retained in each of 



these ca 
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cgorioc, it is necossai*y to look at the nature of tho toi-j and tho effect 
that a plane of symmetry has on then. 

II - 5«i Added i-lass Ter ms 

A body moving through a. real fluid induces a notion in the other-* 
vise stationary fluid because tho fluid must Move aside and then close in 
behind tho body* As a result of this notion the fluid possesses kinetic 
energy that it would not possess if the body were not in notion. The added 
mss toms in tho equations take into account tho energy given to tho fluid 
by the body, 

If tho body Motion is steady, the related fluid notion is also 
steady which requires that the lcinetic energy be constant,. If the kinetic 
energy is constant, no work- is being done on the fluid and therefore the 
added nass terns may be emitted. 

If, however, the body is in accelerated notion, there will be 
work done by the body on tho fluid and it will be necessary to retain at 
least some of the added mass terras. 

Work is accomplished by moving a force through a distance , In 
the case of a submerged body, the distance is the distance the body travels 
and the force is the integral over the surface of the body of the pressures 
exerted by tho body on the fluid. This force, in general, represents a 
system of forces and moments acting on tho body which can be obtained from 
equations (2,3 -- i through 6), when the lcinetic energy is varied. 

The lcinetic energy of tho fluid can be represented as a function 
of the six velocity components (u, v, w, p, q, r), A quadratic form of this 
function as given by Lamb (see page 1?2 ref l) is: 
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2T^. - Au 2 Bv 2 Cvr 2 -r 2A 'v.-r •!* 2B‘tm *!- 20" uv 

+ Pp 2 Qq 2 + Hr 2 •!• 2?‘qr 2Q ' rp + 2R’pq 
+ 2Lup + Pilvq 2rurr 

•!* 2F (vr v v;q) 2G (irp •!• ur) 2H (uq v vp) 

•1- 2F' (vr - vrq) + 20 ‘ (wp - ur) •! 211* (uq - vp) (2,5 - *0 

where the twenty-one coefficients A, B, C etc, are certain constants deter- 
mined by the form and position of the surface relative to the co-ordinate 
axos. 

Letting 

F F* » F F - F* ~ F 2 

II H* = H H - II* b l t 2 

G *!- G* - G„ G - G' 



the last six terms in (2,5 - t) become! 
-1- 2F^vr i G^sp •!- 2H. uq 
+ 2Fpwq Gpur t 2HpVp 



Nine of the coefficients in the above expression may be sot to 
zero if wo take advantage of the assumed xz plane of symmetry. Symmetry 
arguments say that T should not be changed if any of the terras wr, uv, qr, 
pq, up, vq, vr, vrp, ur is replaced respectively by (-v) v, (~v) u, (~r) q, 
(-p) q» ("P) u, (~v) o, C-r) v, (-p) vr, (-r) u* Those terms correspond re- 
spectively with the coefficients A', C‘, P‘, R‘ , L, 1’/, N, G, G‘, 

When T f . is substituted for T in equations (2,3 - 1 thi'ough 6), 
these nine coefficients can be traced to their corresponding terns in equa- 
tions (2*5 - 3), at which time these terns ccn be elii ir . ted fro i the 
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hydro Jynardc oxpi e s sion s » 

Sib: partial derivatives rust be obtained fro;a equation (2*5 ~ 4) 
in order to expand equations (2*3 ~ 1 through 6), these are: 



li£ = Au •»* B'v; v C 'v + Lp + II . q G.r 

V'U 1 2 

— s Bv t A hr *r C'u Kq + H^p 



7>T 



-■ ~ Cw -I- A'v B'u *!- Mr + G p F„q 

/ 7 J . c* 



Bp •!" Q‘r K'q *!• Lu •!- G^ir •! H v 



Dv 

Up 

3P 

rji 

— £ = Qo •!- P 1 r - 1 - R'p + liv + H u F,yr 

12 

d T 

“-f- - Rr •!- P'q -I- Q'p v Nw + B v •<- G u 



(2,5 

(2c5 



(2,5 



(2,5 



(2,5 

(2,5 



*0 



2 

6 ) 

7) 

8) 

9) 

io; 



These derivatives are then substituted in equations (2,3 ~ 1 through 6) to 

obtain : 



'w c 1 


J c 

V 


Lp •; 




•»** G ¥ 
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!■ A'vq 


*. 


B'uq -i 


:• Krq 


+ G 1 pq 


J, 


F 2 qq 






- Ahrr 




C *ur - 


/ 

- Hqr 


~ F rr 
1 




H ? pr 


(2,5 * 


’ 11 ) 


'w + C' 


e 

■u 


4 •; 


:■ F r 


+ h 2 p 










I- B'ur 


•i- 


C * vr * ! 


Lpr 


* 


. 1 . 

» 


G ? rr 






* A'vp> 




B'up • 


• Nrp 


- G,PP 


" 


F 2 q P 


(2,5 * 
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Z s Cu A'v + B'u + Nr + Gp f 2 «| 

Bvp -i- A'wp Cup -!- l v iqp F^rp -!- II^pp 
~ Auq “ B'trq - C'vq - Lpq «• H^qq - G^rq 



(2,5 - 13) 



K » Pp Q'r H*q Lu + 
ltrq •: P'qq •: Q'pq -!- 

- Qqr •• P’r r - P'pr ~ 

•!- Oi rv -i- A'vv •: B'uv 

« B - A'v/.r ■■ C*ir,; - 
vu 

H = Qq -I- P'r v R‘p m •!- 
*!• Ppr -i Q'rr ■■■ P'qr i- 

- Urp " P'qP ~ Q'pp " 

•!- Auu B'trtr •!* C‘vw •!• 
~ Cxro. - A'vu -■ B'uu ~ 

M - Rr * P’q + Kw + F $ * 

*!- Qqp •!* P ’ rp *1- R'pp -!■ 

- Ppq « Q’rq - P'qq - 

Bvu *. u A'vru -1- C ' uu 

- Auv *~ Bhrr •• C'w - 



0 ,v II v 
1 2 

/ 

Nrq F vq 'I G^uq 

Pvr “ II ur ~ F,,rr 
1 2 

lh'7 ’<• G pv *•* F^qv 

Hqir - F rw *• H pw 

H x u •!• Kjr 

Lur G w •: Hvr 
1 2 

/ 

IP-rp ~ F^vp - G^up 

Lp.-r P K^qir ■< G ? nr 

/ 

Kru - G ^ pu - F qu 



Kvp II^up F^irp 

Luo -• G ua - II vo 
1 * 2 

/ 

Mqu *!- F vu + H^pu 
Lpv ** H^civ - G rv 



(2.5 ~ Ik) 



(2.5 » 15) 



(2.5 - 16; 



Sotting the various coefficients in equations (2.5 ~ 11 through 
16; equal to their counterparts in equations (2.5 ~ 3) it is found that? 



2 ? 



A' a X a - X « Y. = - Y « Z. « Z r, K = - K a - M a I! 

vq vrr vr vp v vp vv V ru vu wu 



C' r- X, - - X a Y. « Y a Z = - Z a ~ K a M a I! 

v ur u vu up vq wr vu uu 



w 



G 1 “ X rv, = - Y nn = Z » - * . = K a H 

1 FI PP P vr pv vrr 



~ H . a ~ H 
pn vrq 



G 2 * X i*- = Y rr 



Z a K a 
rq uq 



M a l 'I a N, « - 11 
up nr u rv 



LaX, aK. al-l a H a-H a -l'I 

p u ur p.r uq pv 



M a y. a - K a M, a N a 1) 

q QW v V p qu 



N a z. a K S ~ ft a « N a N 0 
u vrq vrp uu vr 



P‘ B K a - K a ft, a ft a K„ a 11 

qq rr r qu q up 



R* a K, a .. K a ft, a-M a II a }] 

q P 1 * p qp pp qq 

/ / 

1 J ~ ft a X a K 

qu vr 



L - N = Y = 1-1 a II 

PU VU pU 



ft ~ L = Z 



pq 



G, -!• G a 1-i a - ft 
1 2 ru pu 



ft - L a H ~ N 

vp uq 



(2 «5 - 1?) 



H, -1- II a u rr - l'I 
1 2 U P Vc] 



B ~ A = K 



uv 



(2,5 - 18) 
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a « x. = y *= - z 
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vr v vp 
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COlit. 



vhcro equations (2.5 •• 1?) represent the toms vhich arc zero due to 
symmetry. 

Eliminating those terms vhich, duo to tho xz plane of sym otry, 
aro zero, and using the standard nomenclature, as given in reference (6), 
in place of Lamb's notation, the added mss expressions become: 
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The equalities given in equations (2,5 ~ 1 7 and 18) are based 
entirely upon potential theory and do not necessarily hold in tho presence 
of circulation and viscous effects. If circulation acid viscous effects are 
neglected, these equalities provide reasonably good estimates of some of 
the second order coefficients which are not directly amenable to measure- 
ment by conventional towing tank techniques, 

II - 5,2 Second Order forms 

The second order terms which are not added mass terms are : 
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(2.5 - 30) 
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The symmetry plans (xs plt.no) force and moment torus (X, Z, li), 
which involve products of tho syv. ©try plane velocities (u, w, q) with the 
out of piano velocities (v, p, v) rust be zero for the reason that tho cares 
force or moment rust result if v, p, r are replaced by ~v, -»p , -r. This 
doos not hold true with respect to products of xz plane velocities or to 
products of out of plane velocities. 

Another syrsnotry argument that can be made is that a. symmetry 
plane velocity or products of symmetry plane velocities should not cause 
out of the piano forces or moments. For example, there should be no Y 
force resulting from a v: velocity or even a. combination of w and u ve- 
locities , 

After removing the so symmetry terms from equations (2,5 - 2 5 
through 30 ) the second order hydredyname terms left arc: 
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IX •- 5«3 I-’ near Tern s 



The linear 


terms 


in tho hydrodynamic force 


equations are 
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In order to dc terrain© which of those terras should be retained t each terra 
raust, in general be considered on its c;m merits* There is, however, one 
group of terras all of which may be eliminated on the basis of the assumed 
xz piano of symmetry. 

Tho coefficients which, due to symmetry, raust be zero are: 

(1) those that involve derivatives of the symmetry piano forces 
and moment with respect to the out of plane velocities (v, p, 
r), 

(2) the cut of plane force (I) and moments (K, H) which would 
arise from an in tho symmetry plane velocity (u, w, q). 

Those that fall in the first category are of the type which re~ 
quare tho force or moment to stay the same while tho velocity can change 
sxgn. Those of the second category would require a force perpendicular lo 
the piano of symmetry to result from 0 . flow in the pit no of symmetry. 
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Eliminating those terms from the linear terms equations , results 
in tho foil owing '• 
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Discus: 


si. on 


of Terris 






The 


hydrodynamic 


force and moment equations 


that result 


after 



eliminating the symmetry terms will now be looked at in order to further 



reduce tho number of terns in the equations* 

Of first importance is the decision to retain all terns, or a 
form from which they can be derived, that are included in "The Standard 
Equations of Motion for Submarino Simulation" (soo ref 8) that mould apply 
to an unpropolled vehicle. 

In general, the terns which are retained herein but arc not re- 
tained by either NSKDC or 1 -jIT /JL are done so with the idea that not all 
submersiblos possess the near fore and oft and symmetries that modern 



military submarines and the DSRV possess, 

Tho procooduro to bo followed will be to look at tho terns the t 
remain after the symmetry torms have been eliminated ar.d the terms retain- 
ed by other authors have boon set aside. 
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Prior to loolring at the individual feree a;:d no:.- nt ^ual^oiis, 
the effect of the choice of expansion point for the. Tap 3. or seitico should be 
investigated, 

The Host common operating point about vh'j.cb hydrodynamic forces 
and moments aro expanded is some finite forward velocity, than this is 
done, the linear terras given by equations (2.5 ~ 43 through 48) exist. 

These terms include the effect of circulation which does not appear ’in the 
potential theory and the flunk moments which arise from the potential theory. 
Additionally these terms account for the effect that some finite initial 
velocity has upon the a rag terms . 

for the present study of emergency ascent trajectories, the ini- 
tial operating condition is for the vehicle to l -* at rest :ln the fluid. 

Since the force causing the vehicle to ascend is the result of releasing 
ballast from any position on the vehicle, it is as likely for the vehicle 
to start moving astern as it is ahead. l.’ith this sort of operating condi- 
tion, the most reasonable expansion point for the Taylor series is then the 
sero velocity condition. This, however, presents the problem of completely 
eliminating all the linear terms from the series expansion, since they must 
be evaluated at the expansion point, r.oro velocity, for examples 



3v 



= 0 



u ~ u. 



since u G - 0 



Di spite the elimination of the linear terns in th. Trylor scries 
expansion, wo arc dealing with a real fluid and the effect of circualvion 
’.•ill still arise as the vehicle cor: , cnees to novo and the potential ihcory 
strll indicates that the bunk moments exist. The logical place for them to 
bo included is, of course, in the second order terms such as l < , There- 

\\r 

fo-.'c, in the devr Xopm-nl that follows it must 1 •» rcr.mi ucd, ih. t v ]■>.- •<: 5a 
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generally included as a linear effect is nor part of the second order 
effects , 

2 

The drag effects represented by Z r, Z r ' , etc, in the usual 

* i.T \ T V7 

r\ 

case r.iust nor be represented by only tho second order terms Z \j c , etc. 

vn-r 

In the interest of developing as general a set of oquations as 
possible, the linear terms usually appearing in the hydrodynamic equations 
shall be retained, though zero for tliis particular case. 
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The axial force equation without symmetry terms i s : 
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= X.u •: X 0 q 2 X r 2 X vr + X wq v X pr 
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U V XT 

-I- X*u •!- X c q 'l' X uq X \m •!- X vp 

w q uq uvj vp 



*1* X u X xr + X o 
u v q 



( 2 <5 - 49 ) 



where the first tvo lines represent tho terras to be arbitrarily retained. 
Tho last tiro lines contain those terras which require further investigation 
before being retained or rejected. 

Tho linear terra X , will be dropped in favor of tho non-linear 

X which equally well represents tho drag phenomena. This is especially 

u 

true if consideration is taken of tho non-diwonsionalizing parameters in- 

volved. X is non-dimensionalized by dividing by (-g^l u) v/hore U is, in 

general, the velocity of the origin of • the body axes. Therefore, rather 

than leave the dimensions of X^ a function of velocity, we can take a 

further derivative with respect to u and eliminate the velocity dependence! 

This would then give us the alternative non-linear form indicated above , 

The linear terms X and X are assumed to be zero on the basis 

w q 

of experimental results, (See table I of ref 2„) 

The added, mass terms X., X e and X are great3.y dependent upon 

w q uq 

the vehicle shape. If thoro were fore and aft symmetry there would be no 
axial force ro suiting from X, * If, however, the vehicle had a foiiu such 

li- 
as ALVIN, which possess no fore and aft symmetry, there ray well bo forces 



arising flora X, 0 Very much similar arguments can bo said for X, c.r i X 

q 

and therefore they have been i c-tain d iu tho equations. 



uq 



The second order tern is assumed to be zero on the basis of 

experiment;.! results (see table 111 of ref 2), 

Tho one remaining coefficient X is a second order roll trans- 

vp 

verse velocity coupling coefficient which appears to bo essential! y scro, 



II - 'j J\ ,2 la teral F or ce 

The lateral force equation vithout symmetry terms is: 
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whore the first three lines represent the terms to be arbitrarily retained. 

On the basis of the potential theory developed in section II -» 

5.2, Y and Y are of the order of magnitude of X, and X, respectively, 
up ur w u 

Strumpf, in reference (2), retains these terms but notes that there is 

little or no experimental data available for Y while experimental results 

up 

for Y indicate that it is important, 
ur 

Similarly, Y r is also believed to be of importance on the basis 
uv 

of experimental results cited in reference (2). 



II “ 5 « 0 3 No rmal Perce 

The normal force equation without, symmetry terms is: 
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(2c5 ~ 51) 

uhero the first three lines are the arbitrarily retained tonus. 



The added mass terras Z , . Z . Z and Z are retained on the 

u uq 'KJ qq 

possibility of a loss symmetric vehicle’s giving rise to this sort of terra. 

Potential theory estimates Z to be of order X. , which is not negligible, 

uq u ° ° 

Z. and Z to be of order X c , which v.s. s not neglected previously and Z 
u nq V r qq 

of order X, according to potential theory. 

q 

The Z terra is shown to be important on the basis of e: peri- 
uu 

mental results (see ref 2). 



II - 5,4.4 Rol ling Moment 



The rolling moment equation without symmetry terras is i 
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whore the first three lines represent the arbitrarily retained terras. 

Here again there is little or no experimental evidence available 
from which estimates of K^, and K may be. made. Potential theory in- 
dicates that these terms are small and, therefore, three toi c are neglected. 
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3.1 - 5. / K5 Pitching n o : nt 

The pitching Moment equation nit bout symmetry terms 5ss 
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whore the first three linos contain the arbitrarily retained terms. 

Equations (2,5 ~ 18) show that the added mass terra m, is of the 

u 

same siae as X c sand M vrhich have been retained, Similarly ii is of tire 
q qw u q 

same size as Ii. , Experimental evidence cited by Strurapf supports the re- 

vr 

tontion of H , 
uq 

The terra M is, in accordance with equations (2.5 " 18), equal 
uw 

to (X. - Z,), For most vehicles the X, term is considerably D.arp;cr than 
u vr vr 

X, and therefore the term i'l should be retained, 
u uw 

11 - 5«^c 6 Yawing h orrant 

Tlie yawing moment equation without symmetry forms is.* 
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where tho first four lines contain the tews arb'j trarily retained. 



Equations (2,5 “ IS) show that H is equal to Y, •• X , , H " l-I 

vu v Xl up V q 

and N ~ N, . Since Y« is of order of magnitude larger than X«, !•’. is of 
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order of magnitude larger than X, , h is of order Y. and not negligible 
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h , and N end retained by the other authors and therefore, N and 1! 



V Vq 

should also be retained in order to be consistent, 
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II ~ 6 Equations of Motion for Free Asce nt 

The conservation of momentum equations, equations ( 2,3 ~ 32 
tlu'ough 3?), as derived in section II - 3? the gravity forces, equations 
( 2.4 - 2 through 2,4 - 7) from section II - 4 and tho hydrodynamic force 
equations developed in section II - 5 trill now be combined to give the 
final form for the equations of motion for submersible s in six degrees of 



freedom with varying mass and center of mass , 



Several of the terms retained are not generally found in most 
developments because the terms are not experimentally or analytically ob- 
tained by present methods, < 

As it is not the purpose of this paper to evaluate additional 
terms for the equations, these additional terms will be retained but set 
to sere in the computation of ascent trajectories. 



The equations of notion for a vehicle with varying mass rill bo 



presented in the following manners 
represents the rigid body dynamics 



the left hand side of the equation 
the right hand side represents tho 



hydrodynamic forces and moments 
The hydrodynamic term 



acting on the body and causing the motions, 
s presented here have b_;cn non- dirensionalia- 



cd in the usual manner as described in reference (6), 
&*j oral foi'ins of the h^urodyni^ijc tern,') 



Typical non- *dj i • en 



are presented in the nor^ncletiirc * 



For sinpDicdty, the pi’incc have t .cn omitted Iron tho ter 



in the cqu 



pro son led » 



The equations of notion for 



freely ascending vehicle v.-ith \ 



3 able mass are:" 
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(2.6 - 5) 



( 2.6 - 6 ) 



II ~ 7 Su cry 

The axis systems to bo used :in the problem arc discussed a .d t’ : 
transforms lions fron the earth fixed to the body fixed axes are developed. 

A derivation of the dynamical equations for a vehicle vith a vary- 
ing mass and center of mss is then made using a Lagrangian formalism. An 
assumption that the reduction in mass due to releasing ballast can be repre- 
sented as a quasi-steady process is discussed. 

The forces acting on the body are then developed by expanding a 
functional representation of the forces in a Taylor series. The terms in 
the series arc then discussed. 

The equations are then non-dinensionalised foe final presentation. 



.’15 



CHAPTER III 



S0LUT3 OH OF THE EQUATIONS OF nOTlcA 



In the previous chapter a system of equations of motion wei*o de- 
veloped for vhich a solution i ust be found in order to obtain the ascent 
trajectory of a vehicle . It is possible that an analytic solution to this 
set of simultaneous non-linear differential equations could be found, how- 
ever, it is much more practical to assume that a stepuise lineal'* approxima- 
tion to those equations. This approach has been used in many simulation 
studies of vhich the DSRV control system, simulation (so© ref 5) and deter- 
mining the effect of hull shape non-linearities (sea ref 4) are just t>-o. 

The equations as developed in Chapter II conform to the notation 
usod by Strurapf (see ref 2) which does not conform with either the hydro- 
dynamic terms format used by the KIT Instrumentation Laboratory for the 
DSRV control system studies and simulation (see ref 5 ), or the standard 
equations of motion for submarine simulation as used by FSRLC (see ref 8), 
In order that available hydrodynamic coefficients be utilised, the notation 
of Chapter II will be modified in this chapter to conform with that of the 
NSRiC equations. The choice of the NSRDC form over the KIT /IL form was 
made for the reason that it is more likely that vehicle coefficients will 
be obtained from 1.5 ADD than from K1T/IL which must also got its data from 
other sources. The KIT /IL notation does possess the advantage of having 
tho dimensional forms of the coefficients independent of velocity. 

Before a stepwise solution to the equations can be obtained, the 
equations must be put into a form which is a' ocnr.blo to cl is tcs’niou-. 
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form: 



Tho \ chicle equations 



of Motion me generally oicyi csscd in the 



V. 

r ]U , _.U -V -i- A — ^ K 

„ -- * - » (Si X m -■■ V X (V ■:- K X h})) f HW ■: > m ( 3.1 - 1) 

I ~r - U x XW - mBG x (§J u x V) •!• n + n (3d -2) 



dt " ”* "" “ 'dt 

dV „ • * * dvJ _ 



whore ~ " u, w, ~ = p s q, r end I is tho inertial tensor, 
dt dt 

To determine ascent trajectories in a stepwise linear fashion 
on a digital computer it is useful to rearrange the above equations so that 
all the derivatives are on the left side of the equations. The resulting 
equations are of the form: 



dV‘ 
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Tho equations to bo solved for and ~ ■ are then given by th 
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In those equations, unlike equations ( 2 e 6 ~ 1 through ?) fro:.; 

which they come, the linear terns YYY . Z ,Z , M . IT , !i , 11 . K 

v p r u v u \r q v p 

and P arc conbinod vrith the non-linear terms Y , Y ,1 . 7, .7, . i , 

r uv up' ur uu «•••" uu 

u-?’ *uq s \ v ’ %• ^ur* ,jL ^- 5 ’I s done for the reason that these ter. s arc 
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inseparable from a model test data reduction vie. 'point . Tho to: ms Z^ j t 

Z , J-I and ri have been added so as to include all of the hSKDJ cccffi- 
vrw |w| v-vr 

cicnts « The term.' X , X . Z ., 7 , K , it and N have boon included 

hi 1 a mFF j rx EFF EH* El* i' 

in order to allow for the possibility of there being some small of fee tor 
forces on the vehicle. This inclusion in no way affects the legitimacy of 
the equations so long as the force dcos not cause the predominance of one 
velocity component . 

With the equations in this form, a stepwiso linear solution can 
be developed which vrill be capable of being programmed for use on a high 
speed digital computer. 



Ill - 2 Stepwise lincai 1 Sol ution 

A stepwise linear solution is one in which the accelerations are 
taken as constant over a given time interval. The accelerations that exist 
over a time interval are determined from the velocities that existed at the 
end of the previous time interval and the total weight removed from the ve- 
hicle. Thus, the right hand side of equations (3.1 - 3 and 3.1 - h) can bo 
determined, at the beginning of each time step. The inertia terns in 
must also be recomputed before each step since the mass is changing in a 
stepwise linear fashion as described in section II - 3. With the foregoing 
information equation (3»1 ~ 6) can bo solved for the body axes accelerations. 
The velocity of the voliicle in the body axis system may then be 

found by: 

V 



. t 
t f 



(3.2 " 1) 
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where t is tbo time at which t ' u ascent started and t is the present time, 

o 

The integral may, however, be represented as a sen of integrals over each 
step in the stepwise solution. Thus; 
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( 3.2 •• 2 ) 



where t 



h - 1 

and t are respectively the tines at the beginning and end of 



n - 1 

the interval o The solution h„.s been specified as stepwise linear and the 
accelerations are to be obtained as constants for the duration of an inter- 
val, therefore, the accelerations can be removed from the integral leaving 
only the trivial integration of dt from t , to t « Defining a t as 
t equation (3.2 - 2) becomes; 
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( 3.2 - 3 .) 



» 1 

which gives the translational and rotational velocities of the vehicle in 
the body axis system 

j n order to obtain the actual position of the vehiclo relative 
to its inertial starting point the velocities obtained above must bo trans- 
formed into the inertial axis system by use of the inverse of the transfer- 

c 

nation matrixes developod in section II «■ 2*3* Ihs velocities x P y, •* s 

t e 

G, y are obtained from: 
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( 3.2 - l\h 

(3.2 - 5) 

Tho position and orientation in tho inertial axis system it then 
determined by integrating the respective velocities from t to t. Here a- 
gain the solution is stepwise and may bo represented, as a sun of the indi- 
vidual steps t In tliis case, however,, tho integrand x,, , yg, otc. is not a 

h ill 

constant c It is f instead s a linear function of tiue since the acceleration 
from whence it came is a constant* Therefore the equation describing the 
vehicle trajectory and orientation becomes: 



1. 




n - 1 



x. 
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Ut 



( 3.2 - 6 ) 



"il 



where x g ~ (>^, f y g , z ) and - (p, 0,^), 

The computation of x g and ^ a'ee tho final computations of the 
step. The procedure is then repeated until some maximum 'tino is reached 
or some predefined position is reached. 

In order to verify the result obtained by the foregoing pro- 
cedure, end independent trajectory determination must be made, 

111 ~ 3 One Demon s'' onel A s cent Traje ctory 

Now that a three dimensional second order method of determining 
free ascent trajectories has been obtained, some circle on tho results of 
this procedure is in order. The most logical approach is to reduce tho 
problem to one in tho vertical plane only and thereby reducing tho com 
flood ty of the problem. If a derivation of the equations of motion 
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y 



more 



made at tliis point tlio resulting equations would bo no diff exert from tlx^e 
prov.iously obtained except that some ox the coefficients uou.D d now be scro «, 
In addition, the solution of the v..o dimensional equations fox- a ecr .plots 
trajectory would requiro the aid of a high speed digital computer. 

To avoid the necessity of employing a computer, a simplified 
method described by G idling s and Louis in reference ( 9 ) will be used. This 
method, as applied to the problem under consideration, reduces to a one di- 
mensional solution of the equations of motion for a vehicle with varying 
mass. 

The one dimensional solution determines the vertical position, 
velocity and acceleration of the vehicle during its ascentc This inform- 
tion provides a sufficient check of the results of the computer solution of 
the equations of notion developed in Chapter II. A description of the one 
dimensional solution follows 0 

Defining : 



B - the net buoyant force 
b - the time rate of change of B 

t - the sum of all forces in the vertical direction 

g - the acceleration due to gravity 

k - the virtual mass coefficient 

m’ - the mass of the vehicle plus the virtual mass 

W - the instantaneous vehicle weight 

2, 2, 2 - the vertical displacement, velocity and ace 



deration 



respec- 



tively 



rvr 



the crossflow drag and added mass coefficients respectively 



teutons It :s of notion can bo expressed 
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(3-3 •* i) 



The force if' is made up of the voight of the ;jcttiso:incd ballast, 
B, plus or minus the hull drag, D, depending on the direction of notion. 
The force due to the jetty conning of ballast is given by: 



B - \ b (t) dt 



( 3.3 " 2 ) 



lotting the integral bo represented by a series of finite steps the buoy- 
ant force becomes: 
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The mss, m* , of the vehicle can be vrifton as: 



(3.3 » 4) 



(3.3 “ 5) 
(3.3 - 6) 



Substituting these expressions into e.qu' lion (3*3 " ' J ) tho cona- 
tion of notion becomes: 
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Tho technique of step by stop integration similar to that weed in 
III - 2 may now bo applied to this equation in order to obtain the 
nsnsional trajectory. 



cation 
one di~ 



in - 4 



The notation and torus format to be used in 
tion is discussed. The notation used in the equation 
fied to conform, with that of KSR3X and the equations 
ital computer solution. 

A solution method is developed utilizing a 
nique. It accepts the vehicle velocities, position. . 



tho computer Simula-- 
s of motion is modi- 
are rearranged for dig- 



stcqwn.se linear tech- 
orientation and buoy- 



ancy as initial conditions and using equations (3<-l - 6 through 
it computes the accelerations, velocities, and displacements of 
after a time interval At, Before each step the weight, CG, mass 
of inertia are adjusted. 



3.1 - 12 ) 
the vehicle 
and moments 



Finally a simplified one dimens 



onal method is devised for compar- 



ison with the computer results , 
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RESULTS AH; CO./iLV.O C i 

In Chapter II the equations of notion were derived for a fioely 
ascending" vehicle with varying mss and center of ware* A stepviso linear 
solution suitable for programming on a digital computer was dovolopcd in 
Chapter III, and the actual program is presented in Appendix A, An equation 
for a one dimensional accent trajectory van also developed in Chapter II to 
provide a check for the three dimensional, six degree of freedom, ascent 
trajectories program , A program to solve this non-linear differential equa- 
tion is presented in Appendix D. 

In this chpater, .the computer simulations and results of these 
simulations are discussed, conclusions drain and reoemt lend ati ons for future 
work in this area made. 

IV - 1 R esults of Computer Simulations 

Computer simulations of the ascent trajectories for a vehicle 
similar to the Deep Submergence. Rescue Vehicle were conducted using the 
lim 36O computer of the MIT Information Processing Center. These simula- 
tions v r ore conducted primarily for the purpose of debugging the ascent tra- 
jectories program and determining its operating characteristics and second- 
arily to study the motion tendancies of the DSRV. 



IV - 1,1 Program O p era tin g C } • -> c In) -i s ti c s 

The ascent trajectories program in c.ilculc-.t'iug th. vehicle trajec- 
tories depends upon the assumed time incro -’ent for 5 t c aCeu. : _y, Jn goner:.!, 
t " t.n t liner" :.’v a ssu. .or can 1 c class: d as t ■ the r t r o ! ■ ■ , ! < o 5. , It c * 

oil rr^.ht, Rot' that t ; 1 • laid /• cri's.ovy van "'ll ri ' n t "jest it ht' , 



The matter of just tin right incr : lent requires u discussion of the f r.e tors 
influencing tlv situation and, thercfoie, the fo rmor catego* ics trill be 
tieabt vith first. 

Should the tins increment chosen be too large , the forces acting 
on the v'hielc i. ill no l be damped in a naturcl fashion, A vehicle under 
the influence of a steady force till continue to accelerate until this 
force is overcome by another force, and the computer does not see another 
force until the present time increment ends and the forces resulting from 
the present motion can be calculated. If the time increment is too long, 
the notions calculated by the computer become excessive. The eventual re- 
sult is that the computer run is prematurely terminated. The. actual cause 
of termination can be either program caused or machine caused, 

A run is terminated, by the pro.p am if either the time limit c-r 
the depth limit is exceeded, Since reaching the tine limit or reaching the 
depth limbi in a natural fashion do not constitute being premature, they 
are not considered hero. The depth limit can, ho: 'over, be reached in an 

I 

urmatu? i al fashion as is indicated in Appendix C , section C - 3« 

A run is torrdnated by the computer i.f the rmrtcr of ex Ire ; ie 
valuo calculations becomes excessive, 

Terriiination for the latter cause occurs before the normal output 
can be made and, therefore, the only printed output, received i ill be the 
initial condition printout. Termination for the for or cause results in 
all normal output he:b\ printed, 

\ lien the tiue incron^ist assui ad is too s^ell tie computer round 
off on or do. Anutas the actual aalculatjons n d thr output j r* near b^ltss. 
This th n leaves us vith th category of the ti* *> : isr-e . *rt tairy 
"ail id ,hl n . factors i hich o f.r\b! oi a>: effect-* d by th • tv.. ino: • rrl in- 



dud' ballast removal rale, vehicu.D :iv velocities and ace el rata or?, and ve~ 
hiclular stability both static and dyivwic , 

The ballast release rate in pounds per second (real ti»o) can be 
assumed to be fixed for a given vehicle since, in emergency conditions, all 
ballast vill, in general, be released at as high a rate a.s possible* Even 
if this is not done the ballast reD.cc.se rate vrould normally be specified and 
not left to the needs of the simulation. Further, the program in no way 
affects the ballast release rate* 

The velocities and accelerations, on the other hand are directly 
affected by the time increment used. The change in velocity is determined 
by assuming that the acceleration is constant over the time period an d is, 
therefore, equal to the product of the tire increment and the acceleration. 
The acceleration existing during a time increment is dependent upon the 
forces acting on the body during that increment. The hydrodynamic forces 
are, in turn, dependent upon the velocities computed during the previous 
increment. This interdependence of forces, accelerations and velocities 

i 

can be the cause of premature termination if the tine incremv.it is tco 
large. For instance, a large increment vould cause the acceleration to act 
for too long causing the velocity to bo excessive, nhich in turn xrould cause 
the forces in the next increment to bo excessive, e to , , , , 

The static stability of the vehicle is dependent only upon the 
shift of the CG j dative to the CB, and 5s, therfore, a function of the 
ballast release rate only. The dynamic stability is, Forever , dependent 
upon the vehicular velocities c d accelerations as veil i s the change :‘m 
mass caused by the reduction in ui'd, 

The essence of this disc- c.'jo.i 5 th: l, the Li- m Iucxt n" to 
us 1 for vehicle in free ascent i ill b d'.ctmtrd. 1 - .a c Id .lion of t h . 



CO 



ct d accd! 



a lions , 



1; Hast release rate and t> a expect d seed aliens, Th L..11; r t rolea.se 
rate vil’l indicate the orcLr of l.rgnl.tud and tin acceleration:, v ill refin 
tlie increnon l . 

With 'the proper tine incrcmnis chosen trajectories can bo com- 
puted for a variety of initial conditions. 



IV - 1,2 b SKV ] lotion Ten d ancles 

The ascent trajectories generated by the ascent program are not 
those that roe. Id be followed by the D3RV, since the coefficients used in 
the program cone from DSRV model tests with the propeller runnirm , They d.o, 
however, provide an indication of the ascent trajectories and motion ten- 
dancies of that vehicle . Simulations were made with the body axes initially 
coincident with the inertial cozes and with an initial roll angle imposed. 

The results of the simulations with the body axes initially coin- 
cident with the inertial axes clearly indicate the coupling of pitch and 
surge which bring about the trajee Lories of figure IV, 1 , This motion be- 
comes even wore clear in the velocity plots of Figure XV, 2, A ten fancy to- 
ward a negative pitch angle when dumping both trim tanks into the leservoin 
is also demonstrated. 

The simulations involving an initial roll angle indicate the 
snail amount of roll damping and sir- lie stability associated with this ve- 
hicle, They also indicate a slight side force roll coupling as is evidenc- 
ed in the printout in Appendix C, section. C - 2, Th- lack of si bill lg 
also manifested itself by causing the program to terminate pi o maturely when 
using a time increment that proved successful for the s.cro loll cose. Even- 
tually an increment of one tenth th t u r. without rol l promd ruoo -rsful , 
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IV - 2 C orpin s3 one a^d h- 

The computer simulations i si nj the DdkV chr.rao l eristics have dc i- 
onstratsd the v: or th of this pro 0 ram in determining ill'* notions of a vehicle 
as it ascends under the influence of buoyancy alone • Th°y Ivvc shorn that 
a vehicle such as the DSilV is quite sens: '» rave to roll without the stabilis- 
ing effect of a relatively lar 0 e forward velocity. 

The velocities achieved by the vehicle durin w ascent indicate 
that there is a predominant velocity component about v v icli re could expand 
a Taylor series once the initial acceleration phase of the ascent has been 
passed. The most likely t:u ie for a shift, to this sort of forJ ulation is 
twenty seconds after deballastiiv, has been completed. At that tire the ve- 
hicle appears to have reached a terminal condition in which the surge and 
pitch vary botwon essentially constant limits. 

The simulations have indicated that a vehicle vilh near pm,. feet 
jort end starboard symmetry requires no more than too dr a.sfc.ml cqcatlo; r> 
of notion when there are no side forces present, 

i 

It is recoM, landed that future studies L: r.aoo using coefficients 
obtained from ?«odel tests conducted without propellers runnir p . Also, in 
view of the slc rv ascent velocities f it may prove rcrtVuiile to rewrite the 
equations to include the effects of an ascent propulso r vhich r ay prove 
necessary in order to get a more rapid ascent in oth er then c: m regency si t- 
uations . 
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LXrXvi I roil ASCrkT TI.'AJii.C'l nil 



A *> 1 Genera l s 

Tho computer- program to calculate free recent trajectories is 
written for the It ill 3^0 - 65 digital computer of the ITT Tnf01r-2tJ.cn Free* 
casing Center. The program is written in i'OiliAt! TV (see ref 11 end 12 ). 
The program 1 lakes use of the matrix 1 anipulo.tion sub. ‘OU lines contained in 
version III of Die IBM "Scientific Subroutine Package" (ref 13) and the 
mathematical functions contained in the standard “Library Svlieu tines" (ref 
14 ) . 

The system subroutines required are 2 



library Subroutines 5 

SIB - coir, -tos the sine of an angle 
COS - computes the cosine of ca angle 
TAN computes the tangent of an angle 
SQ.UT - computes the square root of m be- . 
ABS - computes the a 1 , solute value of a ir 1 b r 
Scientific Subroutines ; 

Gi-xPlvD - computes the product of t ro r\ ..trie-, s 
IfLHY • computes the inverse of a jr-itrix 



Extensive use is made of the ILb-ilUT fca tm . of i\ t'J.' IV f<v 
input and output. 

The program is divided into seven parts - tho fill! program and six 



c 



( r j 



SubiOUTIkES 



The iiAIK progra • acts as an c: cr.iivc 1 . ading 
tionSf vehicle parafC'lors and program control jnivc/ ■ . oters 
ious subroutines needed to solve the equations of notion 
results of the coni ut a Lions c 



in Initial condi- 
5 calling the vor~ 
and trriting oat the 



The data read in is contained in five )ti\TJLlS r j.3 ; 1. 



ViiiJTS, VJHlCL a.:;d CCJl'ilL, ana one controlled fomst 



?x v. 



put is again in acooida^ce irith the NAIolIJ STS plus a single 



KlUID, 



.si & The cat- 

-L c c nl J C J I 



POSIT, ccntc-dnir^ all the computed results for the entire loot 
Part of the output froia l-JJN jr a set of paehed card: 
tain ell the parameters from IhTftLIST X'uCOI D, This output can 
input data for another run in order to continue the save trajo. 
is desired. A naxhuori of five hundred steps is plloi •_ d in c ,e 



f cried. 

; a hi eh con- 



bo used as 



tor y if it 
run* 



A - i ,2 Subr outi ne C C PF 

C0j 5I’T ? reads in the non-'dircenrirn-,! coeffiei nts and in Hiutoly 
prints thorn for output. The subroutines then computes di. nre ion.fi:- i;.p 
factors, dinensioiialif.es the coefficients for use in the notion equations 
and prints the di v :onsionolij?.od coefficients for out >t. 

The coefficients are read in and printed usi: .. id: it - t' lST 

1 0 2jX 1 i‘ i > c 

Tho dimensionless cceffiei cuts are replaced in storage by the?, 
dimensional fox:. . 



A - 1.3 Sub-- out-5 >r- Ph ’ ST 

BALA3T computes tho changes in ueirh’c 
and moments of inertia due to tho dropping of b 
This subroxxtino is only or Hod Turin- 



mass, center of <"rovilr\ 
Ida at. 

the initial pi- se 



of tho 



trajectory vliru dcbalh asting takes place « 

BALAST calls subroutine /bATRX s.vr.^ - 11 ^ 1 a trio: generated :lr< 

YulATRX varies only during tho dobell as ting phase of tho ascent* 

A maximum of twenty local-’ ens for ballast vreights is allosed. 

There is no printed output from BALA5T, 

A » 1 ,h Subroutin e A HATR X 

AKATRX confutes the six by six matrix of tie c elevation tern coef- 
eficicnts, AH, and then inverts the array for use in the solution of the 
equations of notion. 

A - 1.5 Subroutine HYDRO 

HYDRO computes the hydrodynamic forcing terms end the gravity 
forcing terns and then sums them, plus any effector forces, to get the 
total forcing terms for the equations of motion. 

HYDRO uses the velocities and angles resulting from the previous 
steps in order to compute the hydrodynamic and gravity force components act- 
ing on the vehiclo during the present step. 

There is no printed output from HYDRO . 

A - 1.6 Subroutine TRAJ.OC 

TKAJBC computes tho translational and rotational velocities of tho 
vehicle in both body fixed and earth fixed, axis systems. In eider to do 
tills, TRAJOC also computes tho translational end angular velocity transfer- 
mtion matrices AlhV and TIHV. 

if the pitch angle bccomos ninety ccgi cao the *b an.' for: ution ma- 
trix, A-IliV , rill blow up (mathematically speak h.J. At this t’ , the pro- 
message saying tint this h-; oven- r„d, At the s- .. ti , 

interrupt, Divide ch eld' mo 



gram rill print a 
the computer system will issue a "im 



. :o 



during co. lputation of the Ivo clc cent i in A’L. ’ which involv division by the 
cosine of the pitch angle. At this point the system trill take th.c standard 
corrective action of assigning the value of 10 ; ' : '69 to each of these cle cats. 



A ~ 1.7 Subroutines If SIT '! 

POSITri oojiputes the time ellupsed. and the position of the vohiclc 
relative to the earth fixed origin, The veliicle position includes both 
translation and rotation. POSl'fil also stores all the trajectory gonoi ation 
information in the a may POSIT for the final output sce^uo^oe . 

There is no printed output from POSTTP. 

A - 2 Jn o.it- C’tvl rut 

The variables and constants neediod as input for thus computer 
program arc read into the program by seven road Table 1 lists 

all the input variables and the formats under which they are read, 

A single computer run must be i ;de to a inula to each vehicle. If, 
however, a run tc’ mu. nates before the Vehicle roaches the surface, the 
punched output from the program gives the input necessary to conti hue the 
trajectory on the neat run. 

The program output gives a time history of the simulation. This 
data is stored in an array during execution and printed at the end of the 
run by the main program. During program execution K/u-3IIST / J.KC CM)/ is 
printed at the completion of ei ch step. This is done to aid the user in 
locating the source cf any premature program termination. Premat vo ter- 
mination is generally caused by incorrect input data. 

The dimensions of the storage array rest: -ots fh » prey.- to jjO'J 
til. e intervals, however, the program sto‘. -go >.v .vvbn : --its ■ ■ . , udLh. ih; L 
thus dimension can safely 1 ■> increased to 1000. Thir, coupled with tl e 
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fr om (>> 



ri'U to t r 1 . 



r "I l.r co un - 



ability to continuo a t . ajoctoi » 

3J idled noabor ct iirtci vals during 
For sample. output rco Ap; 






r. trajectory . 
pendix G , 
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TA'-i_- 1 



Ij. - j Ijx *jl '0 ASOhnl hiU'. 1 Ji.i 



Columns Foimut 
2-8 



Symbol 

$irco;!D 



Description 

This namelist mill take more than one cunr-d and 
includes the variables described in Table 2 « 
The list begins x.-ith ylKCCRD end ends v.ith 
#SHD. 



2 - 7 


$FLU1D 


This namelist 
the var 1 able s 


uill take 
described 


one card 
in Table 


2-8 


$VlMTS 


This namelist 
the variables 


v:ill ta.ko 
described 


one card 
in Table 


CO 

6 

Cv> 


$VdI£ECL 


This nr.) olist 


Mill take 


one card 






the variables 


dosenib ad 


in Table 


2-8 


OCOXTRL 


This namelist 
the variables 


v.dll take 
deso:. ibod 


one, card 
in Table 


2 - 8 


$001:13 


This namelist 


L ca.ko 


mere tha; 



includes the variables listed in Tables 2 and 



J o 



1 • 


- 10 


F 10,0 


XI 


11 • 


- 20 


F 10,0 


11 


21 - 


- 30 


F 10.0 


Z'l 



Tha position of tf bo lie. si v.cif 5 t relative to 
the origin of tbo be Jy mis coomli) ntc s< 

(Toet) t 



90 

« “ 



Colui £uS 


Pormt 


Syrifaol 


Fc : r rip tion 


31 - 40 


F 10,0 


P'jASS 


The rate at which the particular bailee u 
weight is removed (pounds per second). 


41 - >15 


I 5 


limi 


The number of tine stops clod re '-'bach a 
particular ballast v,*.i£,iifc is ru.iovcd» The. 
is one card of this type lor each posit! on 
Iron which ballast is re roved « 



Tho folio, -in- tr.W.ob coirt'/fc th 



h/ J.-’iJlo U 0 1..J n.- j^i ^V' 
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TAP. LA 2 
IIAiA'.LISTS 
Paine, Piet /II -‘COLD/ 


PROhiAA USUAL 
VARIABLE KOl-EiroUTUKH 


D3FIKIT10L <?; ULITS 


ISTEP 


The stop number 


DT At 


The time interval used during deball : sing 
(seconds) 


TIU'J t 


The tnrce elapsed since ascent cor.. need (seconds) 


XE x 0 

SJ 


The x component of the distance traveled nod stive 
to the inertial origin ( 3 ecu) 


* 


The y component of the distance traveled relative 
to the inertial origin (feet) 


?.i s. . 

ill 


The z component of the distance traveled relative 
to tbo inertial origin (fret) 


Pin ^ 


The angle of roll (radians) 


THETA G 


The angle of pit -h (ridiens) 



r 


Bo UAL 




VARIAB' 


1IC_ Ji.CJ./'TlJiv' 


DBil 'ill on i luL 


VBARTa 


0 

X 


The vehicle vector of tii lr.t:<c«r:-l v.locih in 
Inertial coordinates (feet per second) 


b’E/BTH 


” o' 

<> 


The vehicle vector of angular velocity in inertial 
coord inc tes (radians per second) 


BVSL 


V 


The vehicle vector of translational velocity in 
tho Body axis sy stoic (feet par second.) 


BROT 


”w 


The vehicle vector of angular* velocity in tho Body 
rods systora (radians per second) 


VJT 


\7 


The vehicle vroipht (pounds) 


XG 


*G 


The x components of tho CG vector (feet) 


l'G 


y G 


The y component of tho CG vector (feet) 

t 


ZG 


Z G 


The z component of the CG vector (feet) 


B 


B 


The vehicle buoyant force (pounds) 


XB 


X B 


The >: component of the CB vector (feet) 


YB 


y B 


The y component of the Cr vector (feet) 


ZB 


f 'B 


The z c O; -pone. it l of the Cf vu-tor (feet) 


1XX 


1 

X\ 


The lass moment of ir.-. rt ’• - .'bout t’x x ads ( py. 


Ti'Y 


J 

yj 


The raoo n. a Bo of in.-, t:» : :• •• ; • y : . is (al.^ 


JZZ 


i 

7.Z 


Tie mass i.. nt cf im! L ->•' ft •. r> v. . 

’T:> 



. :j 



Hi..* : 



. . .l 






fly. 



cc 'J.c - 



lo- List/; ilul)/ 



Ai o jfvav ; c:uj. 



, f / •. *, * .« * 

v Zj-I >• ip'. p <i'l C‘-(/ o > <.> > {.•< 



v / 1 



T ivJv 



ft,. f" dd c.;., ity ( j • r c. d.o feed.) 

fao d.olh of tie voh:Ldc) ;vt it rd. t of the 5 scon 
OVct) 

p a . ) V‘ ■ /vr n.;;/ 

■•• fJXo..ed j'.gc.'. . l ti v - ..■ (;'■. -...'r) 



■br{ 



c.? locftl:a4 f.. i.W .1. b.-.ljv.Gt is i 



3 6£>. r; ccl 



Tho nr 1 L;y- of sic- v rccmiryd i. r - c' ;*loU- cidj<f.laot 



lO 1 



V) PtS 



L:. •, Li.’ / ' • • 1/ 

T(.t 1 • *•; of ’ ' v Jr . 1 r V > ' : 






Voile Jc :u- 



j ' r, ; i 



Nan. lis'c/tt. 







Nan c 1 lise/LC . Vl j v 1 


; 

/ 






iiiOGRAi 

VA21ABLS 


usual 

KU .iLi.u Au**l 


DLFIFITIOF £• It,- 


xs 






X3FF 


Vf 


The X component 


of tl of j. tc. 


fc i’Cto 


(r-^ '-) 


YEFF 


Lff 


The Y co.iipor.ont 


of the ±: A.:;to. 


f urre 0 


(r - .-'..ids) 


ZIFF 


z eff 


The X component 


of tho effector 


,r» . ^ . 

j. one c 


(pounds ) 


KSFF 


1[ EFF 


The K component 


of the effector 


f 01 one 


(p-CUnc! C ) 


KEFF 


l eff 


The in co. ..po-ont 


of the effector 


forces 


(pound •?) 


KFF F 


*Wf 


The N e crap one lit 


of the effector 


• r ‘/ r- 

v^. l> a- >_y 


( P -> . - . r r ) 






Name List/COSFFS 
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X 


X 

■'-j 


Longitudinal io:-. 


‘cv s ec el f?.Or ents .tcv c 


Vohic i 0 






moving ahe ad 








XA 


XA. . 
3-d 


Longi ti. din-:.! f 0 : 

( 


'CO C O <c .1. J. 1L C (j ’ i l-0 


for a 


VL:hjLCl0 



moving as torn 



Y 



Y. . 
13 



Lateral force coefficients for a vehicle movor..^ 
ahead 



YA YA. . 

*3 



Lateral force coefficients for a v< hide moving 
astern 



ii 



Norml force coefficients for a vehicle r.iovJng 
ahead 



ZA 




lor- 'A. force cor ft:' t for a vchid j r 

astern 
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Jiu ) '1 ■ i U>UAL 

ViuO/.l5L3 Al t A '1 lb'. A, JC..^ u t lu 

K K_. . Roll i.Oi .nt x O jX' . c- T. ii./ he .for ; vchloi - j v- 

ij 

inf, ahc : .rl 



IJ- 




H 



I-i. . 



Roll ilo ".‘..'it f o: -oo O' \ .hi ici e to for a vvi 'olo joy- 
i.v.g astern 

Pitch ho. r -'\ :L foj-ct coefficients fee a vehicle 
wovir. ■; alxc d. 



iiA 




Pitch i.ovoiit foi-Ov coefficients for a vehicle 
i.iovlVi 0 astern 



H 



M. . 

i.J 



Yfir iio)..eut fo co coc-^ficit its for a vehicle w- 
ine: chord 



KA 



HA. 



^3 



Yr.vj uO'.cnt fore eo'ffaci-n.ty fox- ?. \ ’'ir."*o r :>v- 



The sttosoidpto i, ;j refer* to the pos; 'bio cc hinc.tj.oi of the 
body velocities and accelerations as described in 'fable 3» 
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Tho stability and control coefficient; vs~d i.; the co. if it or si. - 
uiation or a free ascent trajectory cons fro: 1 a report cn the stability err’ 
control chm-.o-tc:-. id tier of the Deep Sub 'orgerwo he. mo Vehi ole (see ref 15). 
These coefficients irons obtained usinj a no lei ri ih p. speller runrirj, who oh 
violates one of the restrictions originslly imposed on the probl»r . lids* 
however t presents no problen sine;- the eooffioi'V'ts rath, propeller r euri wg 
cm be assigned t' represent sons ficticious vehicle without profiler n:, 
fling# This ficticious vehicle would be si to the i-lh'f but i-cv-ld ef- 

fectively have addition; .3 lifting surface or fin eft. 
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